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Nanopore-based sensing has emerged as a promising candidate for affordable and powerful DNA 
sequencing technologies. Herein, we demonstrate that nanopores can be successfully fabricated in Mg alloys 
via focused electron beam (e-beam) technology. Employing in situ high-resolution transmission electron 
microscopy techniques, we obtained unambiguous evidence that layer-by-layer growth of atomic planes at 
the nanopore periphery occurs when the e-beam is spread out, leading to the shrinkage and eventual 
disappearance of nanopores. The proposed healing process was attributed to the e-beam-induced 
anisotropic diffusion of Mg atoms in the vicinity of nanopore edges. A plausible diffusion mechanism that 
describes the observed phenomena is discussed. Our results constitute the first experimental investigation of 
nanopores in Mg alloys. Direct evidence of the healing process has advanced our fundamental 
understanding of surface science, which is of great practical importance for many technological 
applications, including thin film deposition and surface nanopatterning. 

The last two decades have witnessed an explosion of nanopore technology, which has played a critical role in 
various biological applications, such as the detection and characterization of biomolecules 1 , DNA sequen- 
cing 2 " 4 , and medical diagnostics 5 . Kasianowicz et al 1 pioneered the first in vitro investigations of biomolecule 
transport through single nanopores by passing individual single -stranded RNA and DNA molecules through ot- 
hemolysin protein pores. Recently, attention has shifted to non-biological solid-state nanopores 6 ' 7 and DNA 
origami nanopores 8,9 , which present clear advantages over their biological counterparts 4,5 , including superior 
stability and simple integration into devices. Functionally useful solid-state pores with nanometer dimensions 
have been successfully fabricated in silicon oxides 10 , silicon nitrides 11 , metal oxides 12 , and graphene 13 , via electron- 
beam drilling 11,13 and ion-beam sculpting 14 . 

Because the diameter of a nanopore can be nearly the same length and scale of the translocating species, 
extensive efforts have been dedicated to precisely control the nanopore size, and various methods have been 
utilized 10,12,15,16 . A typical and effective approach is to apply a high-energy electron beam (e-beam), under which 
nanopores present in Si0 2 thin films below a critical size (related to the sample thickness) shrink progressively 
and eventually disappear, resulting in an effective healing process 10 . In contrast, larger nanopores expand under e- 
beam irradiation. Without the e-beam, the nanopore would remain frozen and retain its original size. The healing 
process has been attributed to the surface tension effect of fluidized atoms 10,17 . Using transmission electron 
microscope (TEM) images, which provide real-time feedback during the healing process, allow for the precise 
control of pores with sub -nanometer sizes. 

Although previous investigations on various materials 10 " 13 have propelled the development of nanopore 
technology, to the best of our knowledge, experimental evidence concerning nanopores in metals or metal alloys, 
which exhibit good electrical and mechanical performance, have not yet been reported. Moreover, because solid- 
state nanopores can be employed in biological systems, the compatibility of these materials with biomaterials is 
crucial. Thus, Mg alloys, which possess good electrical conductivity, high mechanical strength and ductility 18 , and 
biocompatibility 19 , were examined in the current study. 

In the present investigation, by applying in situ high-resolution transmission electron microscopy (HRTEM) 
techniques 20,21 , we obtained the first experimental evidence of the healing process of nanopores in Mg. First, the 
nanopores were successfully fabricated in a Mg substrate using a focused e-beam. Subsequently, we obtained 
unambiguous evidence that layer-by-layer growth of atomic planes at the nanopore periphery occurred when the 
e-beam was spread out, leading to the shrinkage and eventual disappearance of nanopores. Thus, the healing 
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Figure 1 | Fabrication of nanopores in Mg. (a-b) Schematic depiction of 
nanopore drilling by a converged e-beam. (c-d) HRTEM images of two 
individual octagonal nanopores with edges enclosed by {0002}, {TlOl}, 
{1100}, and {TlOl} crystal planes. 

process was attributed to the e-beam-induced anisotropic diffusion 
of Mg atoms in the vicinity of nanopore edges. A plausible diffusion 
mechanism that dictates the observed phenomena is discussed in the 
present paper. The atomic-scale observation of morphology altera- 
tions at the nanopore periphery has implications in surface science, 
which is of great practical importance for many technological appli- 
cations, including thin-film deposition and surface nanopatterning. 

Results 

The nanopores were directly drilled in a pure hexagonal-close- 
packed (HCP) Mg substrate (lattice parameter: a = b = 3.21 A, c 
= 5.21 A) by a focused e-beam (10 3 -10 4 A cm" 2 ) inside the TEM, as 
shown in Figs, la and lb. The nanopores were created by removing 
Mg atoms via the knock-on mechanism 22 , which results from local 
exposure to the electron beam (the details of which will be discussed 
later). Depending on the thickness of the sample (10-50 nm), the 
drilling time varied from 5 to 30 s (thicker samples required longer 
drilling times). Interestingly, unlike circular nanopores fabricated in 
isotropic amorphous materials, such as Si0 2 10 and A1 2 0 3 12 , the nano- 
pores in Mg consistently exhibited an octagonal shape enclosed by 
four sets of parallel side- walls (Figs, lb- Id) corresponding to the 
{0002}, {1101}, {1100}, and {IlOl} crystal planes, as observed 
along the [1120] zone axis. Because all of the planes in HCP metals 23 
are typical gliding planes with the lowest surface energies, such octa- 
gonal-shaped nanopores were easier to form due to energy consid- 
erations. The size of the nanopores, which was affected by both the 
converged beam size and drilling time, ranged from 3 to 8 nm (the 
size was calculated by averaging the perpendicular distances of four 
sets of parallel side- walls). Importantly, adjacent atoms around the 
nanopore periphery were well-ordered, implying that oxidation did 
not occur at this stage, in contrast to the results observed in Si nano- 
pores 24 . The obtained nanopores may be applied in molecule ana- 
lysis 25 due to their relatively low surface charge. 

To reduce the intensity to lO^lO 2 A cm" 2 , we defocused the e- 
beam (Fig. 2) and subsequently studied each nanopore via in situ 
HRTEM. Strikingly, the nanopore with an original size of 3.3 nm 
shrank continuously (Figs. 2d and 2e) and finally disappeared 
(Fig. 2f) under e-beam irradiation (see also Supplementary Movie 
SI online), which was indicative of a healing process. The entire 
procedure took only 48 s. When the e-beam was turned off, the 
nanopore retained its shape, suggesting an e-beam-assisted healing 





Figure 2 | Nanopore dynamics under defocused e-beam irradiation. 

(a-c) A schematic illustration showing the healing of a nanopore under 
wide-field electron beam irradiation, (d-f) Time-lapsed images illustrating 
the experimental observations of the healing process of a nanopore with an 
original size of 3.3 nm (Supplementary Movie SI), (g-i) Sequential images 
showing a nanopore with an original size of 6.0 nm under e-beam 
irradiation. The arrow heads indicate the existence of a MgO phase. 

mechanism, which is consistent with the results of previous reports 10 . 
Additionally, the nanopore with a size of 6.0 nm (Fig. 2g), which was 
drilled in the vicinity of the nanopore shown in Figs. 2d and 2e, was 
not able to heal after a long period (701 s) of e-beam illumination. 
Simultaneously, the surface area may have become oxidized, as evi- 
denced by the nucleation of the cubic MgO phase (indicated by arrow 
heads in Fig. 2i). A typical HRTEM image indicating the existence of 
MgO is also presented in Fig. SI. Noticeable expansion of the nano- 
pore was not observed, which may be attributed to oxidation. Such 
size effects suggest that the surface-tension effect 17 should be con- 
sidered in the healing mechanism in Mg. However, the detailed 
relationship between the critical size and sample thickness is beyond 
the scope of the present study and will be studied in future research. 

One distinguished advantage of studying the healing process of 
nanopores in crystalline materials is that it enables the direct visu- 
alization of atomic self-diffusion, which is an outstanding process in 
surface physics 26,27 and is closely related to the intrinsic properties of 
the material. Atomic diffusion also provides one of the fundamental 
mechanisms for the equilibration of crystal surfaces, thin-film depos- 
ition, and surface nanopatterning and is of great value in its own 
right 28 ' 29 . Previous studies on nanopores in amorphous materials 10 , 
which indicated homogeneous contrasts in TEM, may fail to offer a 
detailed mechanism of atom diffusion. Figs. 3a-3f present sequential 
key HRTEM images during the healing procedure of a single nano- 
pore (see also Supplementary Movie S2 online), and Figs. 3g-3l 
provide schematic illustrations related to Figs. 3a-3f, respectively. 
The white dots in Figs. 3a-3f represent actual atomic columns of 
Mg according to the simulated results shown in Figs. S2-S3. Due 
to nanopore symmetry, we only focused on the alteration of the four 
side-walls (corresponding to (0002), (IlOl), (1100), and (IlOl) crys- 
tal planes) throughout the entire process. Compared to the previous 
image, the arrow heads in Figs. 3h-3k refer to newly formed layers. 
For instance, in comparison to Fig. 3g, the growth of the (1101) layer 
is highlighted in Fig. 3h, which altered the surface morphology, as 
indicated by the disappearance of the (1 100) facet (Figs. 3b and 3h). 
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Figure 3 | The healing process of an individual nanopore. (a-f) Sequential HRTEM images showing the alteration of nanopore morphology during 
continuous e-beam irradiation at a constant intensity and focus setting, (g-1) Schematic illustration corresponding to (a-f), respectively. The arrowheads 
point to newly nucleated layers compared to previous images. 



Interestingly, subsequent atom diffusion (Figs. 3c and 3i) led to the 
re-birth of the (1 100) side-wall. Later, the nucleation of (1101) and 
(IlOl) layers again diminished the (1100) facet (Figs. 3d and 3j), 
which re-nucleated, as indicated in Figs. 3e and 3k. Such periodic 
modulation of nanopore surface facets continued until the nanopore 
vanished completely. 

To better illustrate the healing process, Fig. 4a presents the 
observed changes in the size of the nanopore shown in Fig. 3 along 
with the irradiation time. The nanopore width was measured along 
the [0002] (dj, [7704] (d 2 ), [llOOj (d 3 ), and [7704] (d 4 ) directions 
(perpendicular to (0002), (1101), (Il00), and (110l) planes, respect- 
ively, as indicated in Fig. Id), respectively, and the width-time curves 
are represented by black, blue, pink, and red lines. Moreover, Figs. 4b 
and 4c present the width-time curves of typical nanopores with 
original sizes of 4.5 and 3.3 nm, respectively. Initially, the curves 
were characterized by a remarkable stepwise width reduction, which 
was induced by the layer-by-layer growth of crystal planes via surface 
atom diffusion. Each growth event occurred fairly quickly, within 
0.5 s (the limit of our acquisition time). Thus, each step height was 
related to the corresponding inter-planar spacing through multiples 
of integers. For instance, the step height Ad x (Fig. 4a) corresponded 
to a single (0002) inter-planar spacing (2.61 A). The curves consis- 
tently ended with a substantial decrease in the width (indicated by 
arrows in Figs. 4a-4c) (see Supplementary Movies S1-S2 online). 
Such nonlinear shrinking behavior is in contrast with knowledge 
revealed during the shrinkage of nanopores in Si0 2 , where the 
shrinking rate remained constant 1017 . 

Another phenomenon observed in our experiments was that the 
growth of a new (0002) layer typically initiated at a point (indicated 
by arrow heads in Fig. 4a) when the djdi (i = 2,3,4) ratios were 
greater than one (see Supplementary Table SI online). This effect 
was observed in all of the obtained datasets (see Supplementary 
Movies SI -S3 online), except when the original nanopore size was 
small (indicated by a red arrow head in Fig. 4c, where the initial pore 
size was 3.3 nm). Real-time HRTEM images provided in the insets 
(l)-(4) of Figs. 4b and 4c indicated that growth along the [7704], 



[1 100], and [7704] direction was easier to activate than growth along 
the [0002] direction, leading to an increase in the djdi ratios. At 
certain points (d^dj ratios greater than one), the new (0002) layers 
began to nucleate, which often decreased the djdi values to less than 
one. Interestingly, subsequent growth of the (0002) layer did not 
occur until the ratios increased to more than one (i.e., the nucleation 
of a (0002) layer always lagged behind than that of (IlOl), (1100), 
and (IlOl) layers), which was indicative of an anisotropic diffusion 
mechanism and can be explained as follows. Typically, for atom 
diffusion in an HCP crystal, the following two independent jumps 
occur, as shown schematically in Fig. 4d: (A) the jump of atoms to the 
nearest vacancy sites within the (0002) basal plane (indicated by the 
red arrow head), which contributed to the growth of (1101), (1100), 
and (IlOl) layers; (B) the jump of atoms to vacancy sites within the 
adjacent basal plane (indicated by the black arrow head), which 
accounted for the nucleation of a new (0002) layer, as previously 
discussed. Both the theoretical prediction 30 and experimental mea- 
surements 31 indicated that the energy barrier (theoretical: 0.47 eV; 
experimental: 0.57 eV) for jump A is lower than that of jump B 
(theoretical: 0.54 eV; experimental: 0.59 eV). That is, the growth of 
a (0002) layer is more difficult to achieve due to the higher energy 
barrier, which is consistent with the current experimental observa- 
tions. Meanwhile, as the djdi ratios increase, surface tension facil- 
itates the growth of (0002) layers, which reduces the djdi ratios and 
maintains stability. 

Discussion 

The e-beam-induced temperature increase of the sample was esti- 
mated to be less than 1 K (see the Methods section for further 
information) under the current experimental conditions. Com- 
pared to the high melting point of Mg (923 K), such a small temper- 
ature increase is insufficient to explain the fast and spontaneous 
healing of nanopores observed in our experiments. Dynamic displa- 
cements of constituent atoms occurs under electron irradiation at 
high voltages, which promotes atom diffusion in crystalline and 
amorphous phases 22 . The amount of energy transferred by an 
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Figure 4 | Anisotropic surface diffusion procedure, (a) A plot of the pore width versus e-beam irradiation time of the nanopore shown in Fig. 3. 
The nanopore morphology at the moment marked as (b), (c), (d), and (e) is shown in Figs. 3b-3e> respectively. The black arrow heads point to the 
moment before new (0002) layers nucleated, (b-c) Width-time curves of two typical nanopores with original sizes of 4.5 and 3.3 nm. The insets ( l)-(4) 
correspond to the nanopore morphology at the moment marked as (l)-(4), respectively. The width along the [0002], [7704], [TlOO], and [7704] directions 
is represented by black, blue, pink, and red lines, respectively. The arrows point to the observed decrease in width before the nanopore closed, 
(d) A schematic illustration of two independent jumps in an HCP crystal. 



electron to the nucleus of an atom can be estimated according to the 
following equations 22 : 

£ = £ max sin 2 (<9/2) (1) 
£max = £o(1.02 + £ 0 /10 6 )/(465.7A) (2) 

where 0 is the deflected angle of the electron in the field of the atom 
nucleus, E max is the maximum energy transferred, E 0 is the incident 
electron energy (in eV), and A is the atomic mass number. According 
to the equations, the maximum amount of energy (E max ) transferred 
to Mg nuclei by a 200 -kV e-beam is 22 eV, whereas the energy 
required for the displacement of Mg atoms is 9-14 eV 32 . 
Therefore, knock-on damage occurs, which results in nanopore 
formation. In addition, the amount of energy transferred may be 
negligible when 6 is very small. At lower energy beam intensities, 
once the energy barrier along the diffusion path has been overcome, 
atom sputtering is suppressed and nanopore shrinking is driven by 
surface atom diffusion at the vicinity of the nanopore periphery. As a 
result, the healed region (as indicated by arrows in Fig. 2f) with 
brighter contrast is thinner than the nearby areas, as evidenced in 
the simulated HRTEM images (Fig. S4). 



In summary, we demonstrated that nanopores could be fabricated 
in Mg using a converged e-beam. Strikingly, real-time videos directly 
revealed that the healing process of nanopores with sizes of 3-5 nm is 
due to the layer-by-layer growth of lattice planes, which is related to 
the anisotropic self- diffusion of Mg atoms. The diffusion process can 
be initiated by surface tension- driven mass flow, which is activated 
by e-beam irradiation. The ability to fabricate nanopores with con- 
trolled sizes in Mg alloys makes the proposed method a viable 
alternative for potential applications in nanopore technology. 
Furthermore, the healing properties of Mg alloys have extended 
possible applications in irradiated environments, such as in outer 
space. Our work indicates that TEM may be an alternative to other 
techniques, such as scanning tunneling microscopy, in the race 
towards comprehensive investigations of surface science. 

Methods 

The nominal composition of the Mg alloy in the present experiment was 86 wt.% Mg, 
9.0 wt.% Al, 1.0 wt.% Zn, and 4.0 wt.% Sn. First, we made a slice with a thickness and 
width of nearly 500 um and 3 mm, respectively (Fig. S5a), which was mechanically 
thinned to approximately 60 um via mechanical grinding (Fig. S5b). Finally, the 
specimens were ion-milled using a Gatan precision ion polishing system at 5.0 kV. 
During this stage of the process, only the central part of the slice where the nanopores 
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were directly drilled was thinned to the nanometer scale (Fig. S5c). Further details on 
the sample preparation method can be found in our previous paper 33 . The in situ 
HRTEM observations were performed in a JEOL JEM-2010 FEF (UHR) electron 
microscope (Tokyo, Japan) operated at a voltage of 200 kV. 

The temperature increase under irradiation was calculated by Fisher's model 34 : 

AT = ( — - )ln — , where J is the beam current, k is the thermal conductivity, e is 

7TKe d r 0 

the electron charge, b is the sample radius, r 0 is the beam radius, and AE is the total 
energy loss per electron in a sample with a thickness of d. Because energy losses in the 
sample were small compared to the initial energy (200 keV), the term AE/d was equal 
to the stopping power of electrons, dE/dx, which was calculated from the Bethe-Bloch 
equation 35 : 

E(E + mcW 



dE 

"dx : 



27iZp(e 2 /47is 0 ) 2 



2Pmc 2 



(l-/l-|>-' + l>-')ln2+^(l 



where Z is the atomic number of the target element, p is the atomic density, s 0 is the 
dielectric constant, m is the electron rest mass, v is the electron velocity, c is the speed 
of light, E is the electron energy, I e is the average excitation energy of electrons in the 
target, and /? = u/c. 

In the present study, the acceleration voltage was 200 kV; thus, j8 = 0.6946 and 
v = 2.0837 X 10 8 m/s. For Mg, Z = 12, I e = 8.8Z = 105.6, m = 9.3 X 1(T 31 kg, 
p = 4.44 X 10 28 m" 3 (the mass density is 1,770 kgm" 3 ),e = 1.6 X 1CT 19 C,s 0 = 8.85 
X 10~ 12 F/m, melting point = 923 K, and thermal conductivity K = 156 W/mK. The 
maximum temperature increase was estimated to be less than 1 K with 1 = 2.1 nA 36 , 
b = 1.5 mm, and r 0 = 100 nm. 
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